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Abstract: In this paper we deal with the generalized Gamma processes and their com- 
positions. For the compositions of two or more than two generalized Gamma processes 
wc give. wIk^ii possible, the explicit law whereas, in the other c;ases the representations in 
terms of Fox's H-functions are given. We also study the connections between iteration and 
product of random processes by exploiting the properties of the generalized Gamma pro- 
cesses, such a study allows us to obtain some striking result about the compositions of the 
Cauchy processes or fractional Brownian motions. Furthermore, we find out the partial 
differential equations governing the generalized Gamma processes and their compositions. 
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1 Introduction 



A generalized Gamma random variable has a density law given by 

Q{x;c,iJ,,^) = 7 ^^' ^^^^^ , x>0, c,M,7>0 (1.1) 

where c is a scale parameter and n, 7 are the shape parameters. 

The generalized Gamma random variable is very important in the fields of reliability 
analysis and life testing models. It generalizes a number of distributions such as Weibull, 
Raleigh, folded normal, negative exponential. It is well-known that the Gamma random 
variable G has a density law given by the formula (1.1) with 7=1, say Q{x; c, /j., 1). 
In literature, it refers to the Gamma process r{t), t > that is governed by the density 
law 

t — l —X 

q{x,t) = ^^^, x>0,t>0. (1.2) 

The process r(t), t > is a subordinator, a Levy process with increasing paths and 
independent Gamma increments. Its density law is given by the formula (1.1) with the 
shape parameters jj, = t and 7=1. In place of the process T{t), t > we will consider, 
throughout the paper, the processes G-y{t), t > and G^{t), t > which are shghtly 
different from the process T{t), t > discussed previously. The process Gj{t), t > has a 
density law given by the formula (1.1) with the shape parameter c = t. The distribution 
of the process Gj{t), t > is q{x,t] fi,^) = (5(a;;t^/'^,/x, 7) where Q is that in (1.1). 
However, all the processes in question are termed Gamma processes in the sense that 
Gi(l) and r(l) possess a Gamma distribution and Gj{l) and G-y(l) possess a generalized 
Gamma distribution. 

We are interested in the study of the generalized Gamma processes Gj{t), t > 
and G~,{t), t > i] because they generalize a lot of well-known processes as the Brownian 
motion, the Bessel process starting from the origin, the exponential process and so on. 
Furthermore, the compositions of generalized Gamma processes are connected to the 
Bessel functions, in particular, that of the sec;ond kind also known as the Macdonald 
function. It appears truly curious the role that the modified Bessel function Ko{z) plays 
in the product of random processes. 

We show that the distribution q = q{x,t) of the Gamma process Gi{t), t > solves 
the second order partial differential equation with non constant coefficients 

|. = x^5-(M-2)|:<z, x>0,t>0. (1.3) 

Furthermore, for the generalized Gamma process G-y(t), t > 0, we prove that the density 
law, say q — q{x,t), solves the following partial differential equation 



d_ 

Ft 



where 7, /x > 0. From the p.d.e. (1.4) we obtain the partial differential equation governing 
the process G^{t), t > which possesses distribution Q{x;t'^ It suffices to consider 
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that 



d d 

—Q{x;t'^,IJ.,'y) = —Q{x;z,iJ.,-f) 



dt 



(1.5) 



In the present we study the composition of Gamma processes and we provide, when 
possible, an explicit form of their density law, otherwise representations in terms of the 
Fox's H- functions are given. For the composition G^(G^(f)), t > we prove that 



Its density law is given by 



q{x,t;'y,n) 



27 



(f) 



(1.6) 



x>0, i>0, /i>0, 7^0 (1.7) 



where Kq is the modified Bessel function with an imaginary argument. 
We find out a very striking result about the composition G}y{G^^{t)), t > which pos- 
sesses the same distribution of the process G^^{G'^{t)), t > 0. The distribution in question 
reads 



g(x,t;M,7) = 7^^7-^^^p2(^, x>0<>0, 7>0. 



(1.8) 



Such processes generalize, in some sense, the Cauchy process G{t), t > as it appears in 
(1.8). We also show that the following equivalences in distribution hold 



GUGlJt)) 'd- G\{r^)GlJti), i > 0, 7 ^ 



and 



G^^(G2(t)) Gi(t4)G2_^(<5), f > 0, 7 7^ 0. 



(1.9) 



(1.10) 



The last two equalities in distribution show the result mentioned before about the distri- 
bution (1.8). Indeed, the processes (1.9) and (1.10) coincide and therefore they have the 

same distribution. 

We also derive from the distribution (1.8) some interesting result about the Student's 
t-distribution. 

For n independent generalized Gamma processes Gi^. {t), t > 0, j = 1,2, ... ,n we prove 
that 



Gi^ {G^i. . . G!;„ it) . . .)) 'd- H Gl^^ {ti), i > 0, 7i ^ 0, i = 1, 2, . . . , 



(1.11) 

and the distribution can be expressed in terms of H-functions as follows 

, x>0, t>0. (1.12) 



q{x,t) = ^^^Hr'° 



r"(M) 



0,n 
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In the special case where lJ-= \ and 7 = 2, the function (1.8) becomes the Cauchy density 
as we noted above. For the composition of the Cauchy process we prove that 



n 

aC\\C\...\C^{t)\...)\)'^-^C^ ({at)^), t > 0, a > 



(1.13) 



where C^{t), t > 0, j = 1,2, . . . ,n are independent Cauchy processes. 



The process G}^^{G^^{t)), t > has a representation of the density law in terms of 
Fox's H functions which reads 



^1/7172 



(m,0) ; (/x,0) 



, x>0, t>0 



(1.14) 



where 71, 72 > 0. In the case where 71 = 72 we have the distribution (1.7) whereas in the 
special case where 72 = 1 the distribution has the explicit form 



q{x,t;iJ.,-fi) = 2 



71a;' 



.M71-1 



-Kn 2 



a; > 0, t > 0, /i,7i > 0. 



(1.15) 



The iteration of n independent Gamma processes Gi^.{t), t > 0, j = 1,2 ... ,n leads to 
the process 



im = {Gl{. . . G;'+\ (t) . . .)), i > 0, 7,- ^ Vj 

which is related to sonic very appealing results. 

We also show that the following equivalence in distribution holds 

n 

Gi^(G?(...G?(t)...))'=^-nG'^7i(*")' *>0 



(1.16) 



(1.17) 



where I71 1 > and 72 = 73 = • • • = 7n = 1- 

The parameters 71 and /Ui turn out to be the most important parameters and one can 
ascribe such a matter to the fact that Gl^^{t), f > is the guiding process. We remind 

that G'^^(t), t>Q posseses distribution Q(.t; /ii, 71). 
The distribution of the process (1.17) can be expressed as 



q{x,t\ 11) 



7.7Ai-l 



tl/7 



(0, ;^)i=l,2,...,n 



, a; > 0, i > 



(1.18) 



where 71 = 7 > and /ii = > 0. 

We also examine the process B{G\{t)), t> Q where B{t), t > is a standard Brownian 
motion and Gi{t), t > is a Gamma process with distribution Q{x; t, fi, 1). The process 
B{Gi{t)), t > has the distribution 



q{x,t; 11) = 2 — 



, X eR, t>Q, ^l>Q. 



(1.19) 
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Such a function solves the p.d.e. 

A^^q={2f,-3)^q-x^q, xeR\{0},t>0. (1.20) 

where we used the notation q = q{x, t\ /i) for the sake of simplicity. 

It must be noted that compositions involving a fractional Brownian motion BH{t), t > 
and a generalized Gamma process G-y{t), t > turn out to have the same distribution of 
the process B{Gi{t)) for 7 = 2H. In particular, we have 

BH{G2H{t))=- B{\G2H{t)\''")=- B{G,{t)), t>0. (1.21) 

Compositions involving the Gamma subordinator T{t), t > are extensively studied in 
literature (see Madan, Carr, Chung [10] and Madan, Seneta [11] ). 
We study the n— dimensional process 

{B\G^{t)),B'{Gi{t)),...,B^{G^{t))f , t>0 (1.22) 
which is distributed as follows 

^(^'^) = ^^U) ^ ^-t(lNVlj' xeM",t>0 (1.23) 

We give the p.d.e. governing the n— dimensional process in (1.22) which reads 
d 

4— 5= (2At-n)Ag- A(x- Vg), x e M", t > (1.24) 

where A = Yl1=\ '^i is the Laplacian operator, V = Y^^=i is the nabla operator and 
q = qix,t). 

Finally we study a general form of the Bessel process involving the generalized Gamma 
process. As already pointed out the generalized Gamma process G^{t), t > possesses a 
density function which depends on the parameters /i, 7 and c = g{t) (for some positive g). 
This generalization allows to consider a large class of processes for which several results 
are proved to be true. We use the fact that 

i_ 

' n ^ ^ 

E{^7.,mW}''| =G,Mt), t>0. (1.25) 

2 Preliminaries 

We first introduce the Mellin transform of a sufficiently well-behaved function f{x), x € 
(0, 00). The Mellin transform of the function / with parameter 77 G C is defined as 

poo 

M{f{-)}{v)= / x^-'f{x)dx (2.1) 
whenever A4 {/(•)} iv) exists. The inverse Mellin transform is defined as 

•^(^)=9:^/ M{f{x)}{v)x-^drj (2.2) 
^7^* Je-ioo 
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at all points x where / is continuous and for some real 6. 

Let us point out some useful operational rules that will be useful throughout the paper: 



M{fiax)}irj)=a-^M{f}{v), 
M{x''f{x)}{v) = M{f}{v + a), 

M{f{axP)}{v)=^a-"pM{f}('^ 



M 
M 



iv) = (-1) 



{v)=M{f} {v)-M{g} iv) 



r(r7 — n) 



M{f} iv-n). 



(2.3) 
(2.4) 

(2.5) 
(2.6) 
(2.7) 



The formula (2.6) is the well-known Mellin convolution formula which turns out to be 
very useful in the study of the product of random variables. 



Let us introduce the Fox's H-functions. The H-function is a very general function 
defined as follows 



TTm,n 



(ai,ai)i=i,..,p 



where 



nr=im-^/3.)nr=ir(i 



mi) 



j=m+l 



r(i-6, + r?/3,)nLn+ir(%-^«i) 



(2.5 



(2.9) 



and C is a suitable path in the complex plane (for more details, see Mathai and Saxena 
[13]). 

Furthermore, the Mellin transform of a H-function reads 



JO 



?7— 1 TTm,n 



dx = M^;-{v) 



where 



KAni,n 

p,q 



iv) = ^q 

1 Li=' 



nr=i + vpj) nr=i r(i - a, - mt) 



,j=m+i r(i - bj - vl3j) ULn+1 r(«i + mi) ' 

In the sequel we need the following properties of the H-functions 



TTm,n 



{bj,l3j)j=l,..,q 



= C i?"*'" 
P>9 



(fli, cai)i=i^,,^p 

{bj,Cf3j)j=l,..,q 



c> 



(2.10) 



(2.11) 



(2.12) 



{bj,l3])j = l.,..,q 



TT'm,n 



{b3+cPj,Pj)j = l...,q 



c G 



(2.13) 



We also remind the modified Bessel function of the second kind which can be 



presented in several alternative forms: for Jf{j/} > — 5 



argx < TT 



(2.14) 
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(see p. 140 Lebedev [9]) or, for ^ 0, v = ±1, ±2, . . . 

■r, / N IT I-v(x) — Iv(x) , , 

Ku{x) = -^^ argx <7r 2.15 

2 sinz/TT 

(see p. 108 Lebedev [9]) where 

^^(^) = 2^[2) kwik + u + iy N<°°' la^g^K'^ (2-16) 

k—O 

is the Bessel Modified function of the first kind (see p. 108 Lebedev [9]). 
For u >0 there exists the following representation 



i^^ ^y^^ = ^ (a;i^) s^-^e'^'i^ds, p e R, a; > 0, t > 0, C e M (2.17) 

(see p. 370 Gradshteyn, Ryzhik [5]) and the special case where y = 0, p = 2 



POO 

Jo 



Ko{x)= I s~^e~ts^~'^ ds, |arga;|<| (2.18) 



(see p. 119 Lebedev [9]). The Mellin transform of K^, writes 

A^{^.}W = ^r(| + 0r(|-0, ^{n}>u>Q. (2.19) 

We observe that, for v = Q, the Mellin transform (2.19) becomes 

M{K,}{ri)='^^T^Q, 5RW>0. (2.20) 
Finally, in the following it will be useful the approximation 

K^x)^ ^j^, forx^O, y>0 (2.21) 

and 

Ko(x)f^log^, x^O (2.22) 
x 

(see p. 136 Lebedev [9]). 

3 The generalized Gamma process 

In the literature, several authors worked with the Gamma process T(t), t > which is a 
Levy process with increasing paths and cadlag trajectories (a subordinator). Its density 
law is given by 

" ^T(i~' ^>0'*>0- (3-1) 



7 



Hereafter, we deal with the generalized Gamma processes G-y{t), t > and G-y{t), t> Q. 
The first process has a density law (1.1) and the second one is defined by the density 

Q(a:;c,M,7) = 7 ^^p^^^ , x>Q (3.2) 

where 7 > 0, > are the shape parameters and c = g{t) is the scale parameter for some 
positive function g : (0, 00) 1— > (0, 00). The function Q{x; c, fj,, 7) can be rewritten as 

Q{x; c, II, 7) = ^Q7,M (^) 
where, in explicit form 

The function (1.1) is derived as 

Q{x; c', 11, 7) = ^g^,^ ) = *5(a;; c, /x, 7). (3.5) 



Remark 3.1 The function q{x,t) = Q{x;t, iJ,,j) solves the first order partial differential 
equation 

t^^q{x,t) = -~{xq{x,t)), (0,00), i >0. (3.6) 

Indeed, we have 



g^l{x,t) = --q{x,t) + ^q{x,t) 



and 



\ d 1 X \ Li'y — 1 'T x'^ 

— {xq{x, t)) = q{x, t) { g(x, t) -q{x, t) 

7 ox 7 7[x X t 

= - fiq{x,t) + —q{x,t). 

We note that for a process G-y{t), t > (equivalently denoted by G^^^{t), i > in order 
to indicate a process whose density is Q{x; t, n, 7)) with density law q{x, t) = Q{x; t, /x, 7) 
it is straightforward to show that 

[Gi{t)]^/^'= G^{t), t>0, 7>0 and [G^(i)]T =' Gi(i), i > 0, 7 > 0. 

Furthermore, it can be shown that 

[G^(t)]-i =• G-^{t), t>0, j>0. (3.7) 
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where G_^(t), t > possesses the inverse Gamma density law —Q{x;t, iJ.,—'y), /x > 0, 
7 > 0. It will be useful in the sequel to write down the Mellin transform of the functions 
±Q{x; t, /z, ±7). We have 



M{Q{-;t,n,^)}{v)= ' t'-^, 3?M>1-/Z7, 7>0 (3.6 

=E{G,{t)r-' 

and 



M{-Q{-;t,,i,-j)}{v)= 'r(M) ^W<7M + 1, 7>0 (3.9) 

=E{G.,{t)r-\ 
The Mellin transform of the function (1.1) is given by 

m[q{-; t, M, 7) } iv) =M {Q{x; m, 7)} iv) (3-10) 

= roir^*""'' 3?M>i-/.7, 7>o 

=^{g^(*)} 

and, for completeness 

M [-Q{- ; -7)} (r?) =M {g(x; m, -7)} (»?) (3.11) 

= ru^T^*"" ' ^W<7M + i, 7>o 

=E{G_^{t)} . 

For the Gamma process Gi(t), f > with density law q{x,t) = Q{x;t,iJ,, 1) we can 
state the following theorem 

Theorem 3.1 The function 

g(x,t) = ''|'^^'^y , x>0,t>0 (3.12) 
say q = q{x,t), solves the second order partial differential equation 

|. = .^3-(M-2)|^,, x>0,t>0. (3.13) 
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First proof: The time derivative is given by 

The first derivative w.r. to x is given by 
d 1 



dx^ \ X 

and thus, the second derivative w.r. to the space reads 
dx'^ \ x^ x"^ xt t'^ I 



\ a;^ x< a; V 



2 , , ' ' ^ X- 

Xt ' x \ i ^ i2 J « 



_li-2 d Id 
x dx^ X dt^' 
This concludes the first proof. 

Second Proof: We obtain the same result by performing the Mellin transform of the 

function 

(M- 2) 



First of all, by the formula (2.7), we note that 

&^ 
' dx"^ 

and 



A^{^/(^)} = -('?-l)A^{/} iv-l) 

for a function / G (7^(0, oo) such that M{f} exists. 
By combining the previous results one obtains 



M m (??) + M - 2){7j - 1)M {Q{- ; t, m, 1)} (r, - 1) 

=[by (3.8)] = (ry + /i - 2)(,y - 1) Li!^±ii^t'^-2 

=(^_l) r(^+/-l) ,.-2 



at r(/x) 
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This concludes the second proof. ■ 

The generahzed Gamma process is a ^ —self-similar process. We can easily check that 

Gy{at) =• a^G-y{t), i > 0, a > 0, |7| > 0. (3.14) 
Indeed, by the Mellin transform (3.8), we have that 



r ( + i] 

E {G^{at) ^ = {at) ^ = E {G^{t)Y-^ = E ^^a^ G^{t)^ 



'7-1 



for 7 > 0. 

Moreover, by the formula (3.9), we have that 



,r,-l _ V T 



E{G-^{at)Y-^ = -^(at)^ = E {G-^it)}'"'^ = E {a'^ G-^it)}' \ 

for 7 > 0. 

Thus, the Mellin transforms of the distribution of both members in (3.14) coincide for 
|7| > 0. 

We give now the partial differential equation governing the generalized Gamma process 
G.y{t), t > which possesses the distribution q{x,t) = Q{x\t,^,^). 

Theorem 3.2 The density law of the process G-y{t), t>Q 

1 

2:7A'-lg-^ 

qix,t)= ^^^^^^ , a;>0, t>0, M>0, 7>0 (3.15) 

say q = q{x,t), satisfies the following p.d.e. 



First proof: The Mellin transform of the function (3.15) reads 

M/*(77) = [by (3.8)]=r(^'^+M) 3?M>1-7M- (3-17) 
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We perform the time derivative of the formula (3.17) and we obtain 

/^r;-7-l + 7/z\ /?7-7-l \ n-r-i 
' r \- l^]t 1 



7 V 7 / V 7 

= ;^(^ ~ - 7 - 1 + 7/^)*t(^ - 7) 

The inverse MeUin transform, according to the properties (2.7) and (2.4), yields the 
claimed result. 

Second proof: We first evaluate the derivatives w.r. to x 

ox ox ox V X X t 

= ^ {(7M - l)a;^"'^g - 7xg} 

(1 -7)(/i7- 1) 1 d 7 a; 9 

= x^ « + " ('^^-^^ar'^-i^-^iSx^ 

and 

d (1-7) ^ 1-7 d 

By combining the space derivatives, we obtain 

-X ^^?-(7M-l)^- "'^ = -1^-^19^'^ 



7 a; J 7/7, — 1 7 a;' 



.7 

t t \ X X t 



2 I 



-1 



By observing that 

d jj, x^ 

the proof is completed. I 



Remark 3.2 Consider the Theorem 3.2. For 7 = 1 we obtain the p.d.e. (3.13). For 
7 = 2 and /x = | we have the heat-type equation. Indeed, the function 

Q = Q{x;t,\,'2)=^^, |a;|>0, i>0 
2 VTrt 
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satisfies the p.d.e 

d 1 

We observe that the function Q = Q{x; g{t),^, 7), for some positive function g : (0, cxd) 
(0,00), satisfies the p.d.e. 

1 /^.2-^An_r..„_nA.i-7nl'^5 



Furthermore, for 7 = 2, = | and g{t) = 2t, t > we obtain the p.d.e. governing the 
rf— dimensional Bessel process starting from zero. 

3.1 The n-dimensional Generalized Gamma process 

The multivariate Gamma process 

G^,^{t) = {G}^Jh),...,G;jtn)), h,...,tn>0 (3.19) 
possesses dirtibution 

where x, t G and 0=1 — p, 0</9<l. The simbol (/i)^ denotes the quantity 

(/z)o = l, ^fc^, fc = 1,2,.... 

As a direct check shows all the marginals of the distribution (3.20) are Gamma distribu- 
tions. For p = the distribution (3.20) takes the form 

g(x,t)=7"^^-i"« 



r"(M)nr=ii^' 



n 



7/i — 1 



7- 



<rr(/.) 

n 

For n = 2 the distribution (3.20) can be written as 
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After some calculations, we have 



m t^,k\ r(M + fc) 



where 



F(a,6;c;z)=^^^"M%, SR{c} > JR{6} > 0, |^| < 1 (3.22) 

fe>0 ^^'''^ 



r(c) 
'r(6)r(c-&) 



Jo 



is the hypergeometric function (for more details, see p. 238 [9]). 

Also, for a genaralized Gamma process G'_-y,^(i), t > 0, we can write down 



2 .. 1 



thanks to the property 

F{a, b; c; z) = {1 - zy-''-^F{c -a,c-b;c; z) 
(see A S p. 559 [?]). 

Furthermore, for 7 = 1, it is immediate to show that 

k>0 

=H{lJ. + p)tit2. 

The calculations above can be carried out by rewriting the gamma function as 



r(/i + fc + l) = / u'^+'=e-"du. 
Jo 



(3.23) 



By observing that EGi^i_i{t) = ^t, we obtain the covariance function of a Gamma process 
Gi,^(t), t > given by 

Rl{tl,t2) = IJ,ptlt2, ti,t2>0, M>0, 0<p<l 
which is only positive. 
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Remark 3.3 We point out that a fractional Brownian motion has covariance function 
given by 

which becomes Ri{ti,t2) = cr'^tit2 in the special case where H = 1. 

4 Compositions involving Generalized Gamma pro- 
cesses 

Compositions involving two generalized Gamma processes are now examined. Let us start 
from the composition of two independent Gamma processes Gi^{t), t > 0, j — 1,2 each 
with density (1.1) and c = t. Our purpose is to derive an explicit form of the density of 
the compound process G^(G^(f)), t > 0, that can be also viewed as G}y{\G'^{t'^)p), t > 
where the compositions of the processes G^{t), t > will be discussed later on. 

The process (f)). The density of the process Gl^{G'i^{t)), t > writes 

q{x,t;^, IJ.) =-f'^ / s e~y = > e~y*> ds (4.1) 



X > 0, i > 0, > 0, |7| > 



= [by (2.17)] 

_ 2|7| /'x\ii^ lfx\'< 



where Kq is the modified Bessel function (2.18). The Mellin transform of the function 
(4.1) reads 

M{q{-,t;j,fi)}{v)=E{Gl^{G^^{t))} = ^^^^ ^ t^-\ > I-7/X. (4.2) 

In light of the Mellin convolution formula (2.6) we can argue that 

E {G\{G'^{t))y~' = E (t^)}'"' E {G^(i^ )}'"' (4.3) 

where E iGl^{t2) '> for i = 1, 2 is given by the formula (3.10) and thus, we can write 
down the following equality in distribution 

Gi(G2(t))^^-Gi(a)G2(a), t>0. |7|>0. (4.4) 

Thus, from the relation (4.4) we obtain the following equalities 

Gi(G2(f))^=^-Gi(ti)G2(fi), f>0, 7>0 (4.5) 

and 

Gl^iGl^it)) =■ Gl^iti)Gl^iti), t > 0, 7 > 0. (4.6) 
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The distribution of the processes in (4.6) is given by 



q{x,t;-'y,fj,) = -^(-) i^o ( 2W f- , a; > 0, i > 0, 7 > (4.7) 



which is derived from the formula (4.1). 

The process G}y{G'^^{t)) . Wc consider now the compound process Gl^{G'^^,{t)), t > 
obtained by composing two independent Gamma processes. Such processes are, respec- 
tively, a generalized Gamma process Gj{t), t > and an inverse generalized Gamma 
process G-j{t), t > 0. The density law of the process G^(G?_^(t)), t > writes 



,i^,t;,,,)=,^l (-) ^(-) -^ds (4.5 

poo _ 

/ s-^^^-'e-i^y-ii)~"ds 
Jo 



-7' 



(xt)''^ 
xT^Jij 



-' ix^> + t^r^TH,y ->o>*>o>M,7>o. 

The Mellin transform of the function (4.8) reads 

A^M-,i;M,7)}W=r(^^+Ai)r(^i^+Ai) O < 5R{r/} < 7 (/^ + ^ 

(4.9) 

and then, we can ask ourselves if the Mellin convolution results still hold. It is straight- 
forward to rewrite the transform (4.9) as 

i?{Gi(G%(f))}''"'=i?{Gi(a)}''"'ii;{G%(t^)p\ (4.10) 

Therefore, the following equality in distribution holds 

G}^{Gl^{t))=- G\{ti)G^_^{t"^), t>0. (4.11) 
In the formula (4.10) we used the fact that 

'{G_^(i)}''"' = r(^i^ + M) 3?{r?}<7M + l (4.12) 



is the Mellin transform of the density of an inverse generalized Gamma process G-^{t), 
t > as already mentioned in the formula (3.11). 

The process G}^{&_^{t)), t > in the special case = ^ and 7 = 2 becomes a folded 
Cauchy process |G(f)|, t > and its density law has a Mellin transform given by the 
formula (4.9) which reads 



M{qi-,tr,„2)}{v) = E\C{t)\^ ' = ^^, 0<5RM<2. (4.13) 

sm Tj 2 
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Hence, we can write down 

\C{tr^- Gl,^{G\,^{t)), t>0 (4.14) 
and, in view of the formula (4.11), we have 

\C{t)\ 'd- Gl.{thGl,^.{ti) 'd- ^;'fl^ 'd- t > 0. (4.15) 

The last equalities in distribution are due to the fact that [G^{t)]~^ '= G_-y(|). 
Furthermore, from the chain of equalities (4.15), we can state 



Also, note that the relation (4.16) becomes 
1 

for the symmetry of the density law of the Cauchy process. 



Remark 4.1 From (4.9) formula (4.13) follows immediately. We derive the formula 
(4.13) in a different way. It suffices to evaluate the integral 



x^i-^ , * dx=— I x"-! / e-y^^+'^'^Uy dx (4.18) 



7r(t2 + a;2) 





00 


TT Jo 






"ij 










''-\{ 


'J] 











e "^dy 



2 sin 



r 

2>r]>0 



Moreover, 



/ 



l^-i— ^dx = 2 x"-^—-, ^dx= . ^. (4.19) 

' 7r(t2 + a;2) J„ irit^ + x^) sin (77 f) ^ ' 



We consider now the process G^ {G'^{t)), t > 0. It has distribution given by 



X {xt)^^r^n) Jo 



ds 



=7 



x{x^ + t■yy^' V^in) 
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which coincides with (4.8). 

Wc can obtain the same result by considering the MeUin transform (4.9) and the following 
equality in distribution 

Gl^{&^{t)) =■ GL^(i5)G2(fi), f > 0, 7 > 0. (4.20) 
Indeed, from the transform (3.11) it appears that 

E{Gl^{G'^{t))} =^^^-^E{G'^{t)} (4.21) 
and we obtain back (4.9). 

Hence, given the formulas (4.20) and (4.11) we can conclude that 

Gl^{G^^{t))'d-G'^iGl^it)), t>0, 7>0. (4.22) 
Furthermore, from the formula (4.14), we have 

Gl^{&^{t))=- Gl^{ti)G^^{t"^)=-\C{t)\, t>0. (4.23) 
for 7 = 2 and ^Ji= \. 

The composition of two generalized Gamma processes Gl^ ;^i(*)i t > ^-^id G'^ m2(^)' 
t > where the process Gi^.^^.{t), t > 0, for j = 1,2, possesses distribution given by 
Q{x; t, ■^j,iJ.j) are now examined. 

The distribution of the process G}y^^i{G'^^^^^{t)), f > 0, 7 > can be easily written in an 
explicit form as 

^ ^ '(xT +iT)^i+^2 r(^l)r(/X2) ' P-i-A-^ \ J 

For 7 = 2, /xi = i, /X2 = I and v > the distribution (4.24) becomes 

fM=2 ^rV(\ ^ x>0,t>0. (4.25) 

We point out that /(x, z/s) is the Student's t-distribution on the positive real line. 
From the formulae (3.10) and (3.11) we still have 

G';_.(G2^,.(i))'=^-G;^.(i5)G2_^.(tl), t>0. (4.26) 

Indeed, the Mellin transform of the process in question reads 



=i;{G;.(ii)G%,.(i^)}''"' 
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Remark 4.2 Let us consider the process G}y^{G'^^{t)) , t > where G^^(t), t > and 
{t), t > are two independent generaUzed Gamma processes, we have that 

For the right-hand side we are able to give the following expression for the covariance 

function 

E{G}^^{t'^)G'^^{r^)GlJs'^)Gl{.s^)} 
= E{G}^^{ti)G}^y^)}E{Gl{t'^)Gl{s'^)} 
for 7i ^ 0, 72 7^ 0. 



4.1 Compositions involving n independent Generalized Gamma 

processes 

We study, in a more general setting, the composition of generalized Gamma processes. 
We can state the following result 

Theorem 4.1 For n independent generalized Gamma processes Gi^.{t), t > 0, 7^ ^ 0, 
j = 1,2, . . . ,n we have that 

n 

= n^7,(*")' *>0' l7.l>0Vj (4.27) 



Proof: First of all, we write down the Mellin transform of the function Qj = ±Q{x;t, ji, ±7j ) , 
/i > 0, 7j > 0, j = 1, 2, . . . , n which reads 

M{Q,}{v) = ^^t^-' (4.28) 
where 

if Q j = +Q{x; t, /i, +7,) 
< j = l,2,...,n. (4.29) 

12 j if Qj = -Q{x; t, 11, -7j) 

We have to highlight the presence of the term t'^~^ in the formula (4.28) . Such a term is 
closely related to the Mellin transform of the distribution of iterated processes. It turns 
out to be useful in the Melhn convolution machinery as well as we show below. The 
compound process G}:^_^{G'^^{. . . G'^^{t) . . .)) has a density law which is a composition of 
the density laws Qj, j = 1,2, ... ,n. Such a composition can be written as 

@^=i(3j = / .../ Qi{x,si)Q2{si,S2)dsi...Q„{sn-i,t)ds„-i (4.30) 
Jo Jo 

^ V ■' 

n—1 times 
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and thus, we have 



n—1 times 

oo poo 



■ I Tf^rT^^^ ^Q3{S2,S3)dS2...Qn{Sn-l,t)dSn-l 

Jo IMJ i IMJ 



n—2 times 



n . 

=n 



M Vr(M) 



t— , ^ 0, Vj. 



We can easily ascertain that 



Hence, 



ri-l 



E {g^^ (g^^ (. . . g:;„ (t) . . .)) = £; ^ n ^i,, - ) 

and this concludes the proof. ■ 



The previous Theorem is a consequence of the scaling property of the independent 
process Gl^., 1 < i < n. Since Gj{t) = tGy{l) we have, tanks to the independence of 
the considered processes, 

n 

G}^^{Gl{...G;jt)...))'d-tl[Gi^^{l) 

and the claimed result follows immediately. 

The last result is related to the possibility of writing the distribution of the process 
Gj{t), t > (and its compositions) in terms of Fox's H- functions as shown below. 
For the process G^{t), t > with density law Q{x;t'^ , n,^), we have 



X 




t 





(4.31) 



and, for the compound process G}y^{G^^{t)), t > the distribution can be written as 
follows 



r2(M) 



q{x,t) = Hqo 



(4.32) 
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The distribution (4.1) is a special case of the formula (4.32) with 7^ = 72 = 7. 
Finally, we give the distribution of the n— times iterated process (4.27) that writes 



q{x,t) 



r»(M) 



H, 



0,n 



(4.33) 



The representations (4.31) and (4.32) follow from (4.33) whereas, the last one can be 
verified by considering the following Mellin transform 

E{G}^^{Gl{...G;^{t)...))y~' (4.34) 



'r"(M) 



-M < H, 



I3 



n,0 
0,n 



r"(A^) 

= [by (2.3) and (2.13)] 



I3 



~ ^' :^)j = l,2,...,n 



iv) 



n,0 
0,n 



(A'' h)j=l,2,...,n 



iv). 



The study of the composition of Cauchy processes brings to the following interesting 
result 

Theorem 4.2 For n independent Cauchy processes C{t), t > the following equality in 
distribution holds 



n 

\c\\c\...\c-{t)\...)\)\=-^\c'{t^)\, t>0. 



(4.35) 



First proof: It follows from (4.13). We evaluate, making use of the formula (4.13), the 
Mellin transform of the distribution of the folded compound process 

\c\\cX..\c-{t)\...)\ 



t > 



as follows 

E\C\\C%..\C-{t)\...) 



m-l 



1 



sin (r?f ) 
f 1 



^E\C'{\C%..\C-{t)\...)\)r' (4.36) 
E\C\\CX..\C-{t)\...)\)r' 



sin (r?f) 
1 



[sin (r?f) J 

In addition, for the independence of the processes C^{t), t > 0, j = 1, 2, . . . , n, we have 



E 



'7-1 



sin 



This concludes the first proof. 
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Second proof: For the sake of simplicity we consider only one iteration. The general 
case is a very trivial extension. Making use of the representation (4.23) of the Cauchy 
process the following equality holds 

\C\\C\t)\)\ 'd- (4.37) 

for 7 = 2 and /U = | . 

By applying the Theorem 4.1, we can write the formula (4.37) as follows 

Gi^(G2(G3^(G^(t))))^J-GL^(ti)G2(ti)G'%(ii)G^(ti), t>0 (4.38) 
and thus (by (4.23)) 

Gi^(G2(G3 ^(G^(i)))) *=^-Gi^(G2(ii)) Gl,(G^(t^)), t > (4.39) 
'd-\c\ti)\\C\ti)\, t>0. 

By combining such results we obtain 

\C\\C'm\=-\C\t'^)\\C\i^)\ (4.40) 
and the claimed result is obtained. ■ 

The relation (4.14) is equivalent, in some sense, to the well-known subordination law 
given by 

C{t)=-B{Si{t)), t>0 (4.41) 

where B{t), t > is a standard Brownian motion and Si{t), t > is a 1/2— stable law 
(the first passage time). In terms of the generalized Gamma process G-y{t), t > 0, we can 
rewrite the equivalence (4.41) as follows 

\C{t)\'dG^^^(G_,^J^)), t>0. (4.42) 



Remark 4.3 We point out that the process Gl^{G^^{t)), f, > is very different from 
G\{G'i^{t)), t > 0. Indeed the Mellin transform of the density of G\{G^{t)), t > reads 



which is different from the formula (4.2). Such processes have the same distribution if 
and only if 7 = 1. 

The Cauchy process is a self-similar process in the sense that 

^r^/ X 1 atdx tdz „f^, X 

P{C(a,) < „} = ^^^^^^^ = = P{C«) < -} ^ 
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We can also use the fact that 



E\aC{t)\'^-' = a^-'E\C{t)\^-' = — ^ = E\C{at)\''-\ a > 0. (4.44) 

sin(r/2) 

When we compose Cauchy processes, the composition maintains the self-similarity. In 
particular we have 

a\C\\CX..\C''{t)\...)\)\='\C\\C\...\C^{at)\...)\)l t > 0, a > (4.45) 



n 

J||C^ ((a<)^) |, t>0,a>0. 



It is enough have a look at the Mellin transform 

E\a\C\\C^{. . . |C"(i)| . . .)|)| 1""^ =a'>-^E\C^{\C^{. . . |C7"(t)| . . 

=[by (4.36)] = 



(at)'- 



=E 



sin (jjf) 

r,-l 



, a > 



and thus, the Mellin transform of both members in the formula (4.45) coincides. 
Furthermore, for the symmetry of the Cauchy random variable, we can write down 

n 

aCWC^...\C"{t)\...)\) = l[C^(^{at)^^, t > 0, a > 0. (4.46) 

We give now some other particular case related to the compositions of the independent 
generalized Gamma process Gj{t), t>0 which has density law Q{x;t, fi,^), 7, > 0. 

The process G}y{Gl{t)). The composition of a generalized Gamma process G}y{t),t>0 
and a Gamma process Gi{t), t > has a density law which can be written in an explicit 
form as 

g(a;, /X, 7) =7 / r^ds (4.47) 



Its Mellin transform writes 



r 



E{G\iGiit))Y-'= Y2(^,) SRM>1-7M. (4.48) 

Such a transform can be rewritten making use of the Mellin convolution formula (2.6) as 
E{G^(G?(t))}''-^ = £{G':;(i^)}''"'i;{G^(i^)}'"' (4.49) 
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and thus, we can conclude that 

G\{Gi{t))=- Gl^{ti)G^^{t"^), t>0. (4.50) 

Wc also give the distribution of the process G^^^j (Gf (f)), t > where /ii ^ ^2- It can 
be written as follows 

2 f ''\ ^^^^^ / I — 'y^ 



q{x, t; 7, Ml, M2) = ( -r ) I 2\/ ^ ) , xeD,t>Q. (4.51) 



Prom the formula (2.21) we have 

I^. 7, Ml, M2) - r(Mi)r(M2)i''^'^^^^ " ''^^ 

and D = [0, 00) if 7/i2 > 1 whereas D = (0, 00) if 7/^2 < 1- 

Furthermore, due to the symmetry of the modified Bessel function {K_i, = K^,) we have 
the same result in either case 112 > Mi or Mi > M2 (the set D will be determined by 
7Mi ^ !)• 

The equality in distribution G^,^,(G?,^^(t)) =' ^^.^^^(i^/^) G2^^^(iV2), t > holds as 
well. 

Remark 4.4 We can easily ascertain that 

GliiGlli'^t))=' \B\ti)B^{t-^)\ t>0. (4.52) 

Just calculating both Mellin transforms. We remind that G2,i(2i) = |-B(i)|. 
Furthermore, if the subordinator is an inverse Gamma process, then we have (4.42). 

Remark 4.5 We also note that G\{G^^{t)) G\{\G^{r)\'^), t > 0. From the property 
|G^(f)|T =■ Gi{t), t>Owe have G\{\G^t'^)\'^) =' Gi(Gf(r)), t > 0. 

For the sake of completeness wc also give a representation in terms of Fox's functions 
of the density laws that do not have an explicit form. 



The process G}^^{G'^^{t)). The density law of the process G^^(G^^(i)), t>0 writes 

^71-1 POO 

g(x,t;M,7i,72)= ^,,j^y^ s'^-'^-^-^e-^-^ds, x > 0, t > (4.53) 
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where we supposed /xi = /i2 = M > 0, 7i > 0, 72 > 0. 
The Mellin transform of the function q = q{x, t; 11, 71, 72) reads 

,,,,,, r(^ + M)rfe + M) ^^ 



= [by (2.11)] 
=[by (2.3)] 



(4.54) 



(a* ^j!^j)) (m 7172 '7172^ 



(^) 



^1/7172 



Xt 



X 

i 1/71 72 



-J/J I 'yj )l (a* ■ 



(a^,o) 



7172 ' 7172 ' 



(m,o) 

1 



71 ' 71 ^''^ 7172 ' 7172 ' 



iv) 



iv) 



=[by (2.13)] 

tV7i72 I X ^'^ 



f 1/71 72 



(/z,0); (m,0) 
(a*) (M) 7vy2 ^ 



and thus 

a(a^,t) 



1 

^^1/71 72 



H. 



2,2 



.T 

^1/71 72 



(a*,0); (a*,0) 
(M) (a*! 



a; > 0, i > 0, 



(4.55) 



or 



2.il/7i72r2(;,)^0:° 



^1/7172 



(a*) )i (a*! ) 



7l72 ■ 



x>0, t>0, (4.56) 



with Ai,7i,72 > 0. 

In the special case where 71 = 1 and 72 = 7 we have the process G'}(G^(<)), t > that is 
the composition of two independent Gamma processes Gl{t), t > and G^(i), t > 0. 
The density law of the process G'}(G^(t)), t > can be written as 



q{x,t;^) = 



7a;'' 



1 r 



eM7-A»-l 



e « * 



ds, x>0, t>0, /(i,7 > 0. 



(4.57) 



Such a function has the Mellin transform 

Hence, the function (4.57) can be expressed in terms of Fox's H-functions by performing 
the inverse Mellin transform of the formula (4.58). We have 



A^{9(-,i;7)}(^) 



q{x,t\i) = 



1 

r2(/x) 27ri 

1 0-2,0 

^XtV7^2,2 



— ioo 



r(77 + ;u - l)r ( — ^ + // ) tyx-'idrj 



(/i,0); (/x,0) 

(a^,i); (a^4) 



> 0, f > 



(4.59) 
(4.60) 
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where ^ is a real number such that 9 > 1 — n. 

The function (4.60) can be obtained by the formula (4.55). 

By combining a generalized Gamma process G}y{t), t > and an inverse generalized 
Gamma process G'^^{t), t > we obtain the process G^(G?_^(t)), t > with distribution 
given by 

q{x,t; iJ,,-f) 



At7-1 poo 



i^r2(/i) J, 

The Mellin transform of the density (4.61) reads 



/ s-''-T''-^e-Ve-7Tirfs, a; > 0,t > 0, /i,7 > 0. (4.61) 
Jo 



M{q{- ,t;iJ,-f)}{r]) 



r('^ + .)r(,-^) 



-/X7 < ?7 < At7^ + 1. (4.62) 



For 7=1 and n= 1/2, we have 



2 2 TT sm(f[??-|]) 



2<^<2 



(4.63) 



which is the Mellin transform of the distribution of the process |C'(;^)p. This follows 

from the equality (4.42) by noting that [G^{t)X< '= Gi(t). 

We can give a representation of the function (4.61) in terms of H-functions as follows 



_ 7,-1 
t ^ 



= [by (2.3) and (2.13)] 



t ^ 



r2(^) 



xt~ 



iv) 



iv) 



r2(/i) X 



(i-m4) 



and thus 

q{x,t; iJ,,-f) = 
or equivalently 

q{x,t; /i,7) = 

with > 0, 7 > 0. 

We dwell now on the composition 

7S(t) = G:;^(G^,(...G!;+\(i)...)), i>0 



a; > 0,i > 



(4.64) 



t-'' 0-1,1 



xt-' 



(1-M,^); (m,0); (m,0) 
(m, :^); -; 



x>0,t>0 (4.65) 



(4.66) 



26 



where Gi^. (t), t > for j = 1, 2, . . . , n + 1 are independent Gamma processes with density 

Q{x;t, jjL,^j). The most important tool in the study of such a process, once again, is the 
MeUin transform. For the density of the process /^(i), t > 0, say q{x, t), with = jj, > 0, 
7j > 0, J = 1, 2, . . . , n we have 

The transform (4.67) contains fundamental informations about the structure of the mo- 
ments and the density of the process Icit), < > 0, in particular we can state the following 
theorem 

Theorem 4.3 Let us consider the process /^"^(i) = G^^(G^^(. . . G!;^(t) . . .)), t > 
where Gi^.{t), t > 0, j = 1,2, are independent generalized Gamma processes. If 

7i = 7 > and 72 = 73 = • • • = 7n = 1? then the following equality in distribution holds 

n 

Gl^{Gl{...{G^{t))...)y^\[Gl^{t'^), t>0. (4.68) 



Proof: Fixing 72 = 73 = . . . = 7„ = 1 in the formula (4.67), we obtain the Mellin 
transform of the density law of the process lQ~^{t), t > which reads 



1-1 



M{q}iv)={ ^r(^) ' t^} ■ (4.69) 

We immediately recognize the structure of the product of independent random variables. 
The processes Gi^^ {t), t > 0, j = 1,2, . . . ,n are independent and thus 

( n ''"^ n n r f 2=1 + u") 

E I n G^.. (^) I = n ^ {G'.. wr' = n ^%r^*^ • 

The last formula permits us to write down 

E {Gl^^iGli. . . (G?(i)) . . .))}"-' = e\1[ Gi^^it^) ' 



and this concludes the proof. 



In a similar way as in the previous proof, we can show that the process (4.68) is 
71— self-similar in the sense that 

n 

{G^.(G?(. . . {GUm . ■ .))} =■ n Gi,,{{at)i), a > 0, 71 > 0. (4.70) 
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Indeed, we have 



E 



-a-^E{G\^{Gl{...{G-,{t))...))] 



1=1 j ^ i ''71 
=a 71 ^ 

=[by (4.69)] 




=E lllGi^^iiat)^) 





a > 0. 



Furthermore, we observe that 



G^,(Gl(Gf (. . . (G?(i)) . . .))) 'd- G,, n ' * > 




(4.71) 



where we recall that Gj^{t) ~ Q{x;t, iJ,,ji). 

For the process Gl^^{G^^^^{. . . G'^Jt) . . .)), t > where G^^ (t), t > for j = 1, 2, . . . n 
are independent Gamma processes with distribution Q{x;t, ^ 0, jjLj > 0, j = 

1, 2, . . . , n the following equality in distribution holds 



{Gi^(G^^(. . . G^Jt) . . 'd- G}(G^^(. . . G-Jt) . . .)), i > 0, 71 0. (4.72) 



It suffices to consider the process Z{t) = G^^Xit)), t>0 with density law given by 



for some process X{t), t > with distribution fx{x,t). 

We have aheady seen that [G^(t)]T Gi(i) for 7 > and [G^{t)]-'^ =' G_^(t) for 7 > 0. 

i d 

By combining such results we obtain [G-)-(t)]^ =' Gi{t) for 7 G K \ {0}. The density of 
the process Gi{t) is Q{x;t, iJ.,1) and then the process [Z(i)]''', t > has a density law 
which reads 



The formula (4.74) is the density law of the process Gi{X{t)) = [Z{t)Y'', t > for any 
process X{t), t > Q. 

It is also straightforward to ascertain that 




(4.73) 




(4.74) 



{G,{G\{Gl{. . . (Gi'(t)) . . 0))}^ =-Gz n Gi(i") , * > 




n 



= G7(i^)nGi(t^), t>Q 
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for 7,/3 e M \ {0}. 



It must be noted that the shape parameter 71 turns out to be very important and one 
can ascribe such a matter to the fact that (t), t > is the guiding process. 



Remark 4.6 We have aheady pointed out that G}^{G'^{t)) =' Gi(|G2(r)|^) and thus 



G^(G^(...G"(i)...)) =-G^(|G^(...|G!;(f 



0, t>o 



'^•Gi(G^(...G?(t^)...)), t>0 



7 



= GMnGj(i^) , t>0 



Vi=2 



i d 

where we used the fact that |G-y(t)|''' =' Gi{t) and the Theorem 4.3. 

We give now an interesting representation in terms of H- functions of the distribution 
of the process Gi{t), t > and its compositions. We readily have that 



X 




t 


(0,1) _ 



(4.75) 



where Q{x; t, fM, 1) is the distribution of the process Gi(t), t > 0. 

We are also able to write down the distribution of the compound process G}(Gi(t)), t > 
as 



t^r2(/x) 



H, 



2,0 
0,2 



(0,1) ; (0,1) 



(4.76) 



In general, for the n— times iterated process, we have 



q{x,t) 



~ TTn,0 

^M— 1 



(0,l)i; (0,1)2; ... 
X > 0, 00 



(0,1). 



(4.77) 







where Gq'^{z) is the Meijer's G-function (for the definition, see Mathai and Saxena [13]). 
Furthermore, the distribution of the composition (4.68) can be written in terms of H- 
functions as follows 



q{x,t) 



■ Tjn.O 



tl/7 



(0,i)i; (0,1)2; ...; (0,i). 



, a; > 0, i > 0. (4.78) 
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It suffices to evaluate the Mellin transform of the function (4.78). First of all, consider 



(0, ;^)i=l,2,...,n 



(J7 + 7M-1) (4.79) 



by the property (2.4) and, 



M{q{-,t)}{'n) = 



r«(M) 



MiH, 



n,0 
0,n 



(0, ^)i=l,2,...,n 



(»? + 7M - 1) 



by the property (2.3). We obtain the claimed result by observing that 

' ?7 + 7,u - 1 



MiH, 



n,0 
O.n 



(0, ^)i=l,2,...,n 



(ry + 7M-l) = nr 



7 



(4.80) 
(4.81) 

5 7n 



and the formula (4.80) coincides with the Mellin transform (4.67) with 72 = 73 = 
1 and 7i = 7. 

The representation (4.77) follows from (4.78). 

The Gamma process Gi{t), t > possesses a density that is an infinitely divisible law 
and its Laplace transform reads 



(4.82) 



Wc exploit now the infinitely divisibility of the distribution of the process Gi{t), t > 0. 
It is useful to remind that for a Gamma process Gi{t), t > and a generalized Gamma 
process G^,{t), t > the following hold 



[Gi{t)]'^~' G^{t) and [G^(f)]^ =' Gi(t), t > 0. 



(4.83) 



The sum of n independent Gamma processes still has Gamma distribution being the 
Gamma density an infinitely divisible law. In light of the properties (4.83) we can consider 
a sum involving n independent generalized Gamma processes Gl^{t), t > i = 1, 2, . . . , n 
and we are able to write down the following equality 



E H.At)r = J2G\jt) Gi,„^(i), t > 



(4.84) 



and ^ > 0, 7i ^ for I = 1, 2, . . . , n. 

The properties (4.83) also allow us to write the following relation 



1/7 



Xt), t > 



(4.85) 



. i=l 



and thus we can generalize the Bessel process. The density law of the process (4.85) solves 
the p.d.e. 



dt 



d 



7^ 1^ dx 



d_ 
dx 



Q - {'ynn - 1) 



dx 



x>0, t>0,nGn (4.86) 
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as we shown in the Theorem 3.2. 

A straightforward check is the case where n processes G2i{t), t > are involved. We 
obtain 



n 



j2mt)? = G2,ii2t), t>o 



(4.87) 



which is the Bessel process starting from zero. 

5 Compositions involving Generalized Gamma pro- 
cess and Brownian motion 

In this section we will introduce and study some iterated processes. Several authors 
studied the compositions of processes, the most popular one is probably the iterated 
Brownian niotion which has been introduced by Burdzy in 1993 [1] and then studied 
thoroughly by Burdzy [2], Khoshnevisian and Lewis [6] and other authors. The study 
of the iterated Brownian motion has been motivated by the analysis of diffusions in 
cracks (see Chudnovsky and Kunin [3], Kunin and Gorclik [8]). The iterated folded 
Brownian motion |i?^(|i?^(t)|)|, t > 0, where B^{t), t > 0, j = 1,2 are independent 
Brownian motions, has the same distribution of the process Gl i (G^ i (2i)), t > and the 

distribution can be rewritten in terms of H- functions by looking at the formula (4.55). 

The iterated folded fractional Brownian motion \Bjj_^{\B'}j^^{t)\)\, t > where B;^ (t), 
t > for j = 1,2, are independent fractional Brownian motions, possesses the same 
distribution of the process |S(G'_^ i (2^^^^))!, t > and can be seen as the composition 

Gl i(G\^ i(2t^^2)), t > involving two independent generalized Gamma processes. 

Such a distribution can be written either in terms of H-functions as in the formula (4.55) 

or in the form 



Several authors (see for instance Mathai and Haubold [12]) have thoroughly studied the 
function (5.1) and its representation in terms of iJ— functions. The function (5.1) is the 
standard thermonuclear function in the Maxwell-Boltzmann case, in the theory of nuclear 
reactions. It appears very important in physics and astrophysic;s when we consider the 
probability for a thermonuclear reaction to occur in the solar fusion plasma. 
Generally speaking, the results shown in the previous section about the generalized 
Gamma processes and their compositions can be useful to study a number of well-known 
processes. 

In the literature, there are a lot of papers devoted to the study of the variance Gamma 
process where a random time change is given by a Gamma process T{t), t > with 
distribution given by the formula (3.1) (see Madam, Carr and Chang [10], Madan and 
Seneta [11], Kozubowski, Meerschaert and Podorski [7]). The variace Gamma process is 
also termed Laplace motion since the increments follow the Laplace distribution. The 
variance Gamma process B{T{t)), t > has Laplace transform 




(5.1) 



£{g}(A) = 



1 



(5.2) 



(1 + A)* 
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where q, once again, is that in (3.1). It is well-known that the following representation is 
possible 



B{V{t))'^-V^{t)+V2{t), i>0 

where Fi, r2 are two independent Gamma subordinators with different parameters. 
Here we study the compositions involving the Gamma process Gi{t), t > and in par- 
ticular we begin with the compositions where it only appears as random time. For the 
process B{Gi{t)), t > where B{t), t > is a, standard Brownian motion and Gi{t), 
t > is a Gamma process with distribution Q{x;t, iJ.,1), /x > 0, we give the following 
result 

Theorem 5.1 The process B{Gi{t)), t > has a distribution given by 

q{x,t; fj) = I Q{s;t,iJ., l)ds (5.3) 



/2ns 

which solves the following p.d.e. 
d Li &^ 1 

dt''=2d^''--4d^^^'"'^ (^-'^ 
- ~ q, a: e M, t > 0, /X e 0, - U -,00 



where q = q{x,t; ii). 

First proof: From the property (2.21) of the modified Bessel function we can observe 
that 

and 

lim g(a;,t; ;^) = oo. 

Consider now the functions 

e~ 

p{x, s) = and Q{s,t) = Q{s;t, fi,l). 

V27rs 



We evaluate the time derivative 

dt 



-q =[ by (3.13)] = p(.x',.s) |,s— - (m - 2)-| Q(s,t)d. 

/■oo Q2 roo Q 

J P{x,s)s-^Qds-{ii-2) j p{x,s) — Q{s,t)ds. 
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By performing an integration by parts one obtains 

f°° d d 
-J ^{p{x,s)s} —Q{s,t)ds 

roc g 

-(/"-2) J p{x,s) — Q{s,t)ds 



—q =p{x,s)s — Q{s,t) 



=p{x,s)s — Q{s,t) 



8 = 



-(A*-l)y p{x,s) — Q{s,t)ds 



where the coefRcients of the integral pQ'ds have been summed. We integrate by parts 
once again, 



^q=p{x,s)s^Q{s,t) 



s=0 



d 

- sQ{s,t)—p{x,s) 



s=0 



d_ 
ds 



{p{x, s)) Q{s,t)ds 



+ y s-g^{p{x,s))Q{s,t)ds- {n-l) J p{x,s) — Q{s,t)ds. 



Also we observe that 



d { 

(m-I) j p(x,s) — Q{s,t)ds={^,-l)lp{x,s)Q{s,t) 

roo Q 

-j g^{p{x,s))Q{s,t)d 



and 



d 

p{x,s)s—Q{s,t) 



s=0 



=p(x,s){(M-l)-^}Q(s,t) 



3=0 



=(a*- lMa;,s)<3(s,i) 



s=0 



because 



Pix,s)—Q{s,t) 



e 2s I s^e 



s=0 



= 0, in>o). 



s=0 



Moreover, 



-sQ{s,t) —p{x,s) 



gfi Ig 2t 



s=0 



dx^ ' ' tf'T{n) 



= 0, (m>0). 



s=0 
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where p{x,0) = S{x) is the dirac delta function. 
Thus, 

di^'^Jf *^(^^(^'*))*3(s,i)c?s + Aiy —{p{x,s))Q{s,t)ds 

^4 /-oo ,, a2 



'Adx^ J 2 p{x,s)Q{s,t)ds. 



In the last calculations we used the fact that the function p{x, s) satisfies the heat equation 

d 1 52 

Furthermore, 

1 53 e-t^ , , M 92 

= ~4dx-^Jo ^V^s^^'''^'^'^2dx-^''^''^'^ 

1 53 , 92 

In order to obtain the formula (5.4), further calculations are needed. In particular 

— ixqix,t)) = ^3—^+x—^^g{x,s) 
and thus 

This concludes the first proof. 

Second proof: Let us define the Fourier transform of a function / as follows 

Him = I e-^^-f{x)dx = (e-*^^/) (5.5) 
Jm. 

when J^{f}{(3) exists. 
Consider now the operator 



W) = {(2,-3)^-x|,}/, fecr 



where C|"(M) is the space of the rapidly decreasing smooth functions. 
Evaluate the Fourier transforms 

= - 0^Hfm 
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and 



=3/3^ (e-'3-,/)+/33^(e-/5x,/) 
=3/32^{/}(/3) + /33^^{/}(/3). 



Hence, we have 



where /(/3) = 

The Fourier transform of the function (5.3) reads 

and the time derivative is given by 

For the operator (5.4) we have 

13"^ 1 
= - M Y9t(/3) + ^g,(/3) 

= -/^$at(/3)fl- ^ 



2 (l + if) 



2 

which coincides with (5.6). The proof is completed. 

Third proof: Evaluate now, the Mellin transform of the function (5.3) 

2^ /ryx /ry— 1 \ 



A^{g(.,t;M)}(r?) = ^r(|)r(^L + ^j_l_ = ^^(^), ^{r^}>Q. (5.7) 
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The time derivative of the formula (5.7) reads 



= ^(7? - 1)(77 - 2)(7? + 2/x - 3)*i(r? - 2) 

= \v{v - 1)(^ - 2)f ,(,y _ 2) + ^(r; - l)(r; - 2)(2/i - 3)f - 2) 
1 f 5^ 1 2u - 3 f 9^ 1 

The last calculations can be readily verified by exploiting the properties of the Mellin 
transform (2.7). This concludes the third proof. ■ 

The moments of the process B{Gi{t)), t > can be evaluated from the Mellin trans- 
form (5.7) as follows 



E{B{Giit))f' = 2^t{2k + l) = ^Ti^k+-jr{k + ii)^, kGN. (5.8) 

The process Bz{Gry{t)), t > is now examined. In such a process, a fractional 
Brownian motion Bx{t), t > and a generalized Gamma process Gj{t), t > are 
involved. The density law of the process B^{G-y{t)), f > is given by 

q{x,t)=j ^=s^T-i —-ds (5.9) 

2 



='tim[t) ^--i^NVil' -ei^\{0K*>0' 

where /x > and 7 G (0, 2). 

The Mellin transform of the function (5.9) reads 



We observe that, for ^ = ^, the density law (5.9) becomes 



2 



qix,t) = ^^^Koi^x\^^j, xGR\{0},t>0 (5.11) 

and the parameter 7 is not relevant. In particular if we choose 7=1, then we get the 
result of the Theorem 5.1 and this result still holds for all 7 s (0,2). 
We also note that 

BH{G2H{t)) =■ B{\G2H{t)r) =■ B{Gi{t)), t>0,0<H<l. 
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Moreover, by (4.50), we can write 

B{G,{t))'d-B^(^^^B'(^^y t>0 (5.12) 

where B{t), t > 0, B^{t), t > and B'^{t), t > are independent Brownian motions and 
Gi{t), t > is a. Gamma process with density law Q{x; |, 1). 

Wo define and study the n— dinicnsional process 

i;{t)=[B\G^{t)),B%Gi{t)),...,B-{G^{t))f , t>0 (5.13) 

where B^{t), t > 0, j — 1,2, . . . ,n arc independent Brownian motions and Gi{t), f > is 
a Gamma process. The distribution of the process IJ^{t), t > is given by 

^('''^) = ^^0) ^ ^-t(HV')' x.R",t>0 (5.14) 

where j|x|p = X]r=i euclidean norm. 

In the 2— dimensional case, we have 



where ^ ^ 5- For the case /i = ^ we refer to the formula (5.48). 

We provide the following result concerning the process IJj^{t), t > in the n— dimensional 
case 

Theorem 5.2 The distribution (5.14) of the process I^{t), t > solves the following 
partial differential equation 

d 

4— g= (2/i-n)Ag- A(x- Vg), x e K", t > 0, > 0. (5.16) 



Proof: The distribution (5.14) of the n— dimensional process /^(t), i > can be written 

as follows 

At first we point out that 
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where A = 9| is the Laplacian operator. 

We are now able to evaluate the time derivative of the function (5.14) which is given by 



dt 



{p.-2)^Q\ds 



(7(x,t)=[by (3.13)] = ^ p(x,s)|sJJq 

= J p{x,s)s-^Q{s,t)ds- {^1-2) J p{x,t)—Q{s,t)ds 



(5.17) 



where we used the notation 



Q{s,t) 



Let us evaluate the first integral of the time derivative (5.17). By performing two inte- 
grations by parts we have 



d 

Ii ^ p{x,s)s—Q{s,t) 



s=0 



d 



ds 



d 



p(x, s)—Q{s,t)ds — / s-Ap{-x., s)—Q{s,t)ds 



ds 



and 



Ji =p{yi,s)s — Q{s,t) 



s=0 



s=0 



p(x, s)Q{s,t) 

1 f°° I r°° d 

+ 2^ J Pi^^s)Q{s,t)ds- -A J sp{x,s) — Q{s,t)ds 

where we used the fact that J^p(x, s) = ^ Ap(x, s) as we mentioned above. 
After further integrations by parts, we have 



d 

Ji =p(x, s)s—Q{s,t) 



-p{-K,s)Q{s,t) 



s=0 



+ -Ag(x,t) 



s=0 



-A{sp{yi,s)Q{s,t) 



=p(x, s)s—Q{s,t) 



s=0 

5=00 



r°° r°° d 

I p{x,s)Q{s,t)ds - / s—p{x,s)Q{s,t)ds 
Jo Jo 



■p{x,s)Q{s,t) 



s=0 



-Asp{x,s)Q{s,t) 



s=0 



s=0 



1 /■°° 1 

+Aq{x, t) + -A / s- Ap(x, s)Q{s, t)ds. 
^ Jo ^ 



We now observe that 



P(x, s)s—Q{s,t) 



s=0 



=(m- l)p(x,s)Q(s,i) 



=(/i- l)p(x, s)Q{s,t) 



s=0 
s—oo 



s=0 



p{x, s)-Q{s,t) 



s=0 
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because 



p{x, s)-Q{s,t) 



e 2. 



s=0 



Furthermore, 

Asp(x, s)(3(s,i) 



A[]5(x,) 



0, (/X > 0). 



s=0 



s=0 



(m>o) 



s=0 



where j;(x, 0) = nr=i ^(a^i)- 
Hence, the integral X\ becomes 

Ii = {jj,- 2)p(x, s)Q{s, t) 



s=0 



1 f°^' 

Aq{x,t) + -A sAp{x,s)Q{s,t)ds. (5.18) 
4 ^0 



In order to evaluate the last integral in (5.18) we point out that 



rt n ^2 



^[]p(x,,s)— (^p(a;„s)) . 



8=1 » j = l 



It can be immediately verified that 



Thus, 



d I Xi e 2I 



d 



dx'f \ ^/2^TS ) dxi \ s \/27rs j dxi \ s 7 " 



^ 11 n 

Ap{x, s) = - - E n P^^i' 



i=i j=i 



p{Xi,s) + Xi—p{Xi,s) 



n 

-E 



p{Xj ,s)+Xi—Y[ P{Xj , s) 



= -^p(x,s)-i^x,^p(x,s). 



s ' dxi 

1=1 
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In light of the last result, we can write down 



Ji=(M-2)p(x,s)g(s,i) 
= {li-2)p{x,s)Q{s,t) 



if "5 
+ Ag(x, t)-jA\ nq{x, t) ^' a^^^'"' 



s=0 
s=oo 



i=l 



s=0 



n 1 
+ Ag(x, t) - - A(7(x, i) - ^ A (x • V(z(x, t)) . 



The second integral in the formula (5.17) can be evaluated as follows 

d 



pOO 

= -(m-2) p(x,t)^Q(s,t)ds 



(/x-2)p(x,s)Q(s,i) 



1 



s=0 



(M-2)-A5(x,i) 



where we integrated by parts and used the relation J^p(x, s) = |Ap(x, s). 
Finally, by combining all the previous results, we have 



and this concludes the proof. 



Remark 5.1 Wc focus now on the Theorem 5.2, in particular on the case where the 
shape parameter ^ = f . From the formulae (4.83) and (4.84) we know that 

n n 

J2\B\t)\''dJ2G\^{2t)'dG,^.{2t), t>0. (5.19) 

Hence, we can write down 

B{G,,^{2t))'d- B{BSQl{t)), t>0 (5.20) 

where BSQ^{t), i > is the n— dimensional squared Bessel process starting from zero. 
Squared Bessel processes are Markov processes, and their transition densities are known 
explicitly. For the ^— dimensional BSQg{t), t > starting from a; > the transition 
density is given by 

g(y,i;x) = -(^) ' e-^/^ f^j, y > 0, t > 0, S > 0. (5.21) 

Iiy{z) is the Bessel modified function of the first kind (sec formula (2.16)). 

The function (5.21) reduces to the Gamma distribution when the starting point is a; = 0, 

and we obtain 

a(y,^;0) = g^|^, y>0,t>0,6>0. (5.22) 
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Now, for = ? in the p.d.e. (5.16), we have a new p.d.e. which writes 



4^9= -A(x- Vg), xeR", t>0 
and, for q = q{x, 2t), we have 



d 



1 



A(x- Vg), xeM", t>0 



2 

The distribution 



(.t)tr(t)""Vx/i;' 



X G 



t > 



which comes from the formula (5.14), satisfies the p.d.e (5.24). 

Finally the p.d.e. (5.24) is the governing equation of the n— dimensional process 

qi2t) = {B\BSQlit),B\BSQ'iit),...,B^{BSQ°M)f , i>0 

which possesses the distribution (5.25). 
For n = 1 we have 

I\{2t) = B{G,^.{2t)) =■ si(|S2(t)|2), t > 
and the governing equation becomes 
d 



dt 



1 

'2dx^ 



d 
ox 



a; e K, i > 0. 



Remark 5.2 Consider the process 

ii^;(i)ii 



E 



t > 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



where B^ (t), t > , j = 1,2, ... ,n are independent Brownian motions starting from zero 
and Gi,^(f), f > is a Gamma process. The process (5.29) is the euclidean norm of the 
process'in (5.13), 7^(i), t>0. 
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We can find the distribution of the process t > by observing that 



= [by (4.83) and (3.14)]'^- 



= [hy (4.68)] =■ 



i.d. 
i.d. 



[by (4.84)] ^2Gi,^(a)Gi,5(a), t>0 

[by (4.68) and (3.14)] 'd' ^d,. (2 Gi,^(t)) 

'^■[by (4.72) and (4.83)] 'd' ^2,5 (2 Gi.^t)), t > 0. 

Therefore, the distribution of the process ||/^(f)||, t > follows from the formula (4.51). 
The process ||/;^(i)|| =' ^j2G\ i^{ti)Gl^^{ti), i > can be also written as 

\\i;m'^\/2G\jGl^), t>0 

'd [hy (4.72) and (4.83)] 'd' x/2G^,^(G?,^(t)), t > 
'^■[by(3.14)]'='-G^,^(2Gf,„(t)), i>0 
'd-Glj2BSQ'M'^)), t>0. 
Hence, it follows that 

\\q{2t)fd- \B{BSQl{t))\, t>0. (5.30) 
The process (5.30) can be written as follows 

||/|(2t)|| =11 {B\\B{t)\^),B^mt)\% • • ■ , B'\\B{t)\')f II, t>0 



\B{\B{t)\' 



t > 



and thus, 



iiB(|B(i)i^)ir=^-|B(iiB(t)r)i, i>o 

where B(t) = {B\t) , B'^ (t) , . . . , B"" {t)f , t > 0. 



(5.31) 



Remark 5.3 Wc consider the composition G^ 1 (G^ ^ (t)), t > with distribution (4.25) 
examined in the previous section, where two independent generalized Gamma processes 
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are involved. The process T{t) = Gl 1(6^2 K{t^)), t > 0, ioi j = 2, has a distribution 
given by 

q{x,t) = 2 



x>0, t>0, iy>0. 



(5.32) 



Furthermore, by keeping in mind the property G^^^{t^) '= G^^^{t), i > 0, 7 we have 
T{t)'^G\.{G\.{t-')), t>0 (5.33) 
=-Gl^{\G\.{t-')\'), t>0 
'd[hy (4.83)] =• Gl ^{G\u{t-^)), t > 
'^•[by (3.14)] 'd- GU{2G\.{2t-')), t > 



i.d. 



B 



B 



t > 0, 



t>0. 



(5.34) 



where the process BSQ°{t), t > is the z/— dimensional squared Bessel process starting 

from the origin. 

Furthermore, from the property 



we are able to write down 

T{t'')=Gli{t-^)G\.{t-^), t>0 
'd-G\At)G\At-^), t>0 



(5.35) 
(5.36) 



i.d 



■[by (3.7)] 



' [by (3.14)] '=^- 



G\ 



t > 



i.d. 



B{t) 



t>0 



where BS^{t), f > is the z/— dimensional Bessel process starting from zero. 
We observe that (by the symmetry of the distributions) 



B 



BSQO (^) 



B 



^.d. B{t) 



BSO (i) ' 



t>0. 



(5.37) 



We point out that, for each z/ > 0, the random variable T{i') possesses the density law 
q{x, v) which is the distribution of a Student's random variable with v > Q degrees of 
freedom. 
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We have already studied the process B{Gi{t)), t > where B{t), f > is a standard 

Brownian motion and G'i(i), t > is a Gamma process with distribution depending on 
the parameter /i > 0. A particular case where /i = ^ is now examined. 

Consider the process Bjj^{\Bjj^{t)\"^), t > where B^^ (t), i > 0, j = 1,2 are indepen- 
dent fractional Brownian motions with covariance function given by 

E{B'H,m'HM)} = \ + l^l'""^ - ' 3 = 1,2 (5.38) 

for < Hi, H2 < 1. 

The distribution of the process B}j^{\B'jj^{t)\"i), t > reads 
where, for the sake of simplicity, we set H2 = H . 

The distribution (5.39) comes from the formula (4.47) with 7 = 2 and n = ^ or the 
formula (5.3) with /i = ^ by considering the symmetry of the distribution. 

We give the following result (see, for more details [4]) 

Theorem 5.3 The process Bjj^{\Bjj^{t)\'^), t > has a density law (5.39) which solves 
the partial differential equation 

', = -H,t^-^-^(2l, + .l^)q (5.40) 



dt \ dx^ dx^ 

for allx^O,t>0 and H2 e (0, 1). 

Proof: Prom the Theorem 5.1 we know that the distribution of the process B{Gi{t)), 
f > is solution to the p.d.e. 

d 52 
A-Q= {2ii-5)—Q-x-^Q, xgR, t>0, ii>0. 

In the special case where fi= ^ we have 

d 92 g3 

4-0 = -2j—Q - xj—Q, x€R,t>0. 
at ax^ ax° 

We also know that B]j^ (|B|^ (*) I ^ ) =' ^(Gi {2t^")) where Gi has a distribution Q{x; 2t'^", 
and thus 



dt'-Jt^' ^dz' 



= -iJi^^-M 2-^ + I g, xeK, t>0. 



Moreover, from (2.21), we have to observe that 
lim q{x,t) = 00. 

X— >0+ 



g2 g3 X 



' dx"^ dx^ J 
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The proof is completed. ■ 

For the process Bjj{\Bjj{t)\Ti) =' B^{\Bjj{t)\^), t > wc prove that 

E{Bl,(|B|W|^)}''"' = £;{s^(ii)}''"'i;{B|(i4)}'"' (5.41) 
and thus we can write down 

Bh{\B%{t)\i)'d- Bl{r^)B%{th t>0. (5.42) 

2 2 

We observe that 

{B\\Bj,{t)\')B\\BUs)\')}'^{B],{t)Bl{t)B],is)Blis)} 

L2222J 

(by (5.42)) and 

E\B],{t)Bl{t)B],{s)Bl{s)] =E\B],{t)B],{s)] E \ B%{t) Bl{s)] 

L2 2 2 2 J L2 2JL2 2 J 

4(I*I" + N"-|*-^I")'- (5-43) 

by the independence of the processes Bjj{t), t > and B]f{t), t > 0. 
The covariance function (5.43) can be only positive. 

We have already seen in the formula (5.28) a special case of the p.d.e. (5.40) where 
the subordination is made with respect to the fractional Brownian motion insteed of the 
Brownian motion. 



In a general setting we can state the following result (see, for more details [4]). 

For the process 

ir\t) = Bh{\BU- ■ ■ \Bm\^ ■••)!*), t > 0, e (0, 1) (5.44) 

with B-'fj, j = 1, 2, . . . , n independent fractional Brownian motions, the following equalities 
in distribution hold 



n 

I^-\t) 'd- B\\B\ . . \B^{t)\^ . . .)n 'd- i > 0, if e (0, 1). 

i=l 

The proof of the first equality follows from the property Bu^t) B(t^^) for each t. 
The proof of the second equality can be given by evaluating the Mellin transform of the 
density of Ip~^{t), t > which reads 



3? .2 





n-1 j 


Jo 




r25r(i)l 


n 

lH{r,-l) 




n 


s.r,-l 



2si 



OtSi . . . aSn-i 



si'" 



i=l 
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Such a transform is in line with the Mellin convolution machinery. 
Also, we remind that 

\BHit)f =■ 1^2,1 =• G,^.{g{t)), t > 

where g{t) = 2t^^ . 
Hence, we have 

I--\t) '^B (g^(G^(. . . G-;^{g{t)) ...)), t>0 

= G2,iKi(G^(...G!^,]^(5W)-..)), t>0 
[by the Theorem 4.3] 

n 



i.d. 



n 



n 

'd-l[Bi{t), t>0. 

The last two steps follow from the property Bnit) B{t^^), for each t> Q. 
Prom the covariance function (5.43), by considering the Theorem 5 we obtain 

{n n ^ 

and 

{n n ^ f n n ^ 

n B% (^) n (^) \ =[by (5-45)] = n % n % (^) [ 

1 /, , 2H I I 2H , , 2ff \ " 

We observe that the n— dimensional process 

{B\\B{t)\% B\\B{t)\^), B-{\B{t)\')f , i > (5.46) 
is a special case of the process (5.13), in particular If (2t), t > 0. Such a process possesses 

2 

a distribution given by the formula (5.14) with li= \ which satisfies the partial differential 
equation 

d 

2 — g= (l-n)Ag- A(x- V?), x e M", i > 0. (5.47) 
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In the 2— dimensional case we have the distribution (5.15) which becomes 



1 1 



e-ll-ll 



where we used the property K_\{^z) = Ki{z) = ^ P' Gradshteyn and 

Ryzhik [5]). 

Therefore, it follows that 



p(x,t) =9(x,2t) 

1 e" 



(a;i,a;2) eM^\{0,0}, t>0. 



is the distribution of the process ||/i(2i)|l = J\B^(\B{t)\^)\^ + \B'^{\B{t)\^)\^ , t > 0. 

2 

We give now some connections between the iterated fractional Brownian motion and 
the composition of the generalized Gamma processes. 
Consider the process 

Pn{t) = B [\{ G'^^ (2i) j , i > (5.49) 

where B{t), t > is a standard Brownian motion and Gi^.{t), t > 0, j = 1,2, . . . ,n are 
independent generalized Gamma processes. The process Pn{t), t > possesses a density 
law given by 

q{x,t)^ — =0"(s;2f,M, 71,72, •••,7n)c^s, xGR,t>0 (5.50) 
Jo v27rs 

where Q"{s; 2t, /i, 71, 72, • • ■ , 7n) is the distribution of the product YVj=i ^4 (2^), t > 0. 
Consider now jj = for j = 1,2, . . . ,n. The Mellin transform of the formula (5.50) is 
given by 

M {q{x, t)} (77) = (I) s^-'Q"{s; 2t, /x, 71, 72, • • • , ln)ds. (5.51) 

The last integral is the Mellin transform of the process (4.66) and is given by the formula 
(4.67) which coincides, for /x = |, with the Mellin transform of the distribution of the 
iterated fractional Brownian motion. Thus, we can write down the following equality in 
distribution 

^[n^'^(2t)j 'd B\,^{\Bl^{...\B-ll^{t)\...)\), t>0. (5.52) 



47 



Furthermore, for fi = 1/2 and 7j = 1 for j = 1,2, ... ,n (or equivalently Hj = 1 for 
j = 1, 2, . . . , n in the process (5.52)) we have 

'd-B\\B^{\...\B-+\t)\\..\)\'), t>0 
where Hn+i = H € {0, 1). 
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